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Abstract: Biological soil crusts (BSCs) are bio-sedimentary associations that play crucial ecological roles in 
arid and semi-arid regions. In the Gurbantungeut Desert of China, more than 27% of the land surface is 
characterized by a predominant cover of lichen-dominated BSCs that contribute to the stability of the 
desert. However, little is known about the major factors that limit the spatial distribution of BSCs at a 
macro scale. In this study, the cover of BSCs was investigated along a precipitation gradient from the 
margins to the center of the Gurbantunggut Desert. Environmental variables including precipitation, soil 
particle size, soil pH, electrical conductivity, soil organic carbon, total salt, total nitrogen, total phosphorus 
and total potassium were analyzed at a macro scale to determine their association with differing 
assemblages of BSCs (cyanobacteria crusts, lichen crusts and moss crusts) using constrained linear 
ordination redundancy analysis (RDA). A model of BSCs distribution correlated with environmental 
variables that dominated the first two axes of the RDA was constructed to clearly demonstrate the 
succession stages of BSCs. The study determined that soil particle size (represented by coarse sand 
content) and precipitation are the most significant drivers influencing the spatial distribution of BSCs at a 
macro scale in the Gurbantunggut Desert. The cover of lichen and moss crusts increased with increasing 
precipitation, while the cover of cyanobacteria crusts decreased with increasing precipitation. The cover 
of lichen and moss crusts was negatively associated with coarse sand content, whereas the cover of 
cyanobacteria crusts was positively correlated with coarse sand content. These findings highlight the need 
for both the availability of soil moisture and a relatively stable of soil matrix, not only for the growth of 
BSCs but more importantly, for the regeneration and rehabilitation of disturbed BSC communities in arid 
and semi-arid lands. Thereby, this study will provide a theory basis to effectively increase soil stability in 
desert regions. 
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1 Introduction 


Biological soil crusts (BSCs) are complex matrices of soil particles, cyanobacteria, algae, lichens, 
mosses, liverworts and microfungi (Belnap, 2001). In undisturbed areas, BSCs can completely 
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cover interplant surfaces and thus constitute as much as 70% of living cover in arid areas (Belnap, 
2003). Despite their unassuming appearance, the diminutive organisms that comprise BSCs, 
surprisingly play many important ecological roles in arid lands. BSCs are known to increase soil 
stability (Eldridge, 1998; Eldridge and Leys, 2003) and nutrient availability (Bowker, 2007), to 
influence soil hydrological cycles in arid lands (Belnap, 2006; Garibotti et al., 2018), to improve 
seedling establishment of grasses and shrubs (Deines et al., 2007; Peter et al., 2016), and to 
provide habitats and food sources for reptiles and invertebrates (Lalley et al., 2006a; Li et al., 
2011). Crust function is tied to crust community types (Pietrasiak et al., 2013) and surface 
morphology, and these variables vary greatly between climatic regions (West, 1990). 

Many studies have examined the relationships between environmental variables and crust 
composition and cover to determine which variables have the greatest impact on crust 
development. In the absence of disturbance, many biotic and abiotic factors have been found to 
contribute both to distribution patterns of BSCs and to differences in species composition. These 
factors include precipitation (Belnap et al., 2007; Biidel et al., 2009), geomorphology (Kidron et 
al., 2010; Pietrasiak et al., 2011; Williams et al., 2013; Bu et al., 2016), soil texture (Lalley et al., 
2006a), soil pH (Rivera-Aguilar et al., 2009), dust deposition (Kidron et al., 2000; Lalley et al., 
2006b), vascular plants (Bowker, 2007; Li et al., 2010), and UV-B radiation (Castenholz and 
Garcia-Pichel, 2000). Most of these studies were concentrated on either hot or cold deserts, such 
as the Colorado Plateau (Pietrasiak et al., 2011), Great Basin Desert (Belnap and Lange, 2003), 
Mojave Desert (Williams et al., 2013), Negev Desert (Kidron et al., 2010), Namib Desert (Biidel 
et al., 2009), Kalahari Desert (Thomas and Dougill, 2007), deserts in Australia (Eldridge, 2003; 
Rogers, 2006) and Zapotitlan drylands (Rivera-Aguilar et al., 2009), but little is known in drier, 
temperate deserts, such as the Gurbantunggut Desert, where BSCs are patchy in distribution. 

The Gurbantunggut Desert of China is remote from the sea and is surrounded by high 
mountains. It is a typical temperate desert with a climate permanently in drought. 
Lichen-dominated BSCs are abundant in the Gurbantunggut Desert. Zhang et al. (2004) described 
the characteristics of BSCs distribution on sand dunes in this desert and found that the species 
composition varied with different successional stages of BSCs. Chen et al. (2007) suggested that, 
at the sand dune scale, both thickness and cover of BSCs are linked to soil physical-chemical 
properties and topsoil textural stability. At a regional scale, lichen and moss crusts are found 
mainly in the southern region of the Gurbantunggut Desert. Crust cover becomes sparse to the 
northern part and is absent in the eastern and western parts of the desert (Zhang et al., 2007). 
However, it is not clear how abiotic factors influencing species composition and BSCs 
distribution in this temperate desert. 

In recent years, BSCs in the Gurbantunggut Desert have undergone significant disturbance as a 
result of intensive human activities, particularly grazing, and road building and construction. 
Locations that have experienced frequent disturbance are characterized by extensive areas of 
loose sand and are largely devoid of vegetation (Wang et al., 2006). BSCs are completely 
destroyed by bulldozers for the use of road building and the construction of oil production 
infrastructure, and the ensuing crust removal is followed by significant changes in the structure of 
the soil surface. Although disturbance intensity appears to be the primary predictor of crust 
development, some relatively undisturbed areas cannot develop crusts. Therefore, a macro-scale 
study can help to explain the mechanism of the formation and distribution patterns of BSCs. This 
data can be utilized to improve strategies for conservation of existing BSCs resources, and 
provide effective suggestions for soil rehabilitation and regeneration of BSCs in disturbed areas. 

Our primary goal for this study was to understand the factors driving the distribution and 
abundance of BSCs in the Gurbantunggut Desert. Specifically, our objective was to assess the 
distribution and abundance of BSCs throughout the Gurbantunggut Desert by studying the 
relationship between biotic and abiotic factors with the development and distribution of crusts. In 
the Gurbantunggut Desert, mean annual precipitation decreases significantly from the outer edge 
inward to the center of the desert and this is accompanied by increased evaporation from the 
boundary to the center of the Jungger Basin in which the desert is located. Frequent sandstorms 
lead to a spatial variation in soil grain size at a macro scale, which significantly affects the spatial 
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differentiation of vascular plants (Qian et al., 2009). Therefore, we hypothesized that precipitation 
and soil grain size would be the major limiting factors for determining the spatial distribution of 
BSCs at a macro scale. 


2 Materials and methods 


2.1 Study area 


The study was carried out in the Gurbantunggut Desert (44°11'-46°20'N, 84°31'—-90°00’E; Fig. 1) 
located in the Jungger Basin, Xinjiang, China. This desert is the largest fixed and semi-fixed 
desert in China with an area of 4.88x10+ km?. The desert landscape is dominated by a mix of 
longitudinal, crescent-shaped and honeycomb sand dunes. Elevation ranges from 312 to 743 m 
a.s.l. The desert is characterized by an arid climate. Mean annual precipitation varies from 70 to 
160 mm, falling predominantly in spring. Mean annual potential evaporation is 2607 mm. Annual 
mean temperature is 6°C—10°C and the maximum summer temperature can exceed 40°C. Wind 
speeds are highest in late spring, with an average of 11.17 m/s. The prevailing wind comes from 
the WNW, NW and N directions. The soil stability of the Gurbantunggut Desert is dependent on 
native vascular plant vegetation and well-developed surface cover of BSCs (Chen et al., 2007). 
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Fig. 1 Locations of the five study sites in the Gurbantunggut Desert. Note: each site comprises 120 quadrats. 


2.2 Experimental design 


The field investigation was started in October 2012. In autumn, BSCs are conspicuous and readily 
identified because most vascular plants (shrubs, subshrubs and perennial herbs) have either 
disappeared or returned to low levels of physiological activity. Five study sites (Fig. 1) 
representing five desert regions were selected to evaluate the spatial distribution of BSCs in 
different regions of the desert, and to determine the environmental variables affecting the 
distribution of BSCs at a macro scale. A nested sampling design was used, i.e., six 10 mx20 m 
plots with well-developed BSCs and spaced 500 m apart were established in each of the five sites. 
It should be noted that these plots were distributed between sand dunes to eliminate the 
interference of spatial heterogeneity. Within each plot, four 5 mx5 m subplots were placed in the 
interdune without causing any disturbances. In each subplot, five 50 cmx50 cm quadrats (2.5 
cmx2.5 cm grid; 400 points per quadrat) positioned at 3 m intervals were selected to determine 
the cover of BSCs, giving a total of 120 quadrats at each site (6 plotsx4 subplotsx5 quadrats). 
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Prior to taking measurements, the soil surface was sprayed with water to make the BSC 
components more readily visible. BSC cover was estimated as the proportion of the 400 points 
per quadrat occupied by the visible components on the crusts (lichen, moss and cyanobacteria), as 
described by Li et al. (2010). The cover of BSCs was then categorized as lichen crust cover, moss 
crust cover and cyanobacteria crust cover. Cyanobacteria often occur in the earliest successional 
stages of BSCs, followed by lichen, and finally, moss (Bowker, 2007). Moreover, in order to 
discuss the influence of environmental variables on the succession of BSCs, we defined four 
dominant BSC types by their relative cover: cyanobacteria crusts, cyanobacteria and lichen crusts 
(cyanobacteria dominant with lichen present), lichen and cyanobacteria crusts (lichen dominant 
with cyanobacteria present), and lichen and moss crusts (both lichen and moss present). 


2.3 Soil sampling and analysis 


Three soil samples were collected from a depth of 0—10 cm in each subplot and then taken back to 
the laboratory where each soil sample was air-dried, crushed and passed through a 2-mm sieve. 
Soil pH was determined in a 1:5 soil:water suspension by a calibrated pH meter (PHS-4, Jiangsu 
Manufactory of Electrical Analysis Instruments, Jiangyin, China). Electrical conductivity (EC) 
was measured by DDS-11A. Soil organic carbon (SOC) was analyzed by the K2Cr207 method 
(GB9834-88); total nitrogen by the CuSO4-Se powder diffusion method (GB7848-87); total 
phosphorus by the NaOH Melting-Mo Te Sc Colorimetry method (GB7852-87); and total 
potassium by the NaOH Melting-Flaming luminosity method (GB7854-87). Total salt was 
measured using methods described by the Nanjing Institute of Soil Research, Chinese Academy 
of Sciences (1980). All these parameters were selected for their potential to explain the 
distribution of BSCs (Bowker et al., 2006; Li et al., 2011). Soil particle size was determined using 
standard dry-sieving techniques. Resultant material was then allocated into four categories: coarse 
sand (0.500—2.000 mm), sand (0.250—0.500 mm), fine sand (0.0625—0.250 mm), and silt (<0.0625 
mm). Levels of calcium carbonate (CaCO3) were determined by titration (Morse et al., 2007). 


2.4 Climate data collection 


Precipitation data from 2005 to 2011 were obtained for the five study sites. Precipitation data for 
site 1 was from the Fukang National Field Scientific Observation and Research Station for Desert 
Ecosystems (44°17'N, 87°56'E; 475 m a.s.l.). Data for sites 2 and 3 were obtained from the field 
meteorological stations (Watchdog1000, SPECTRUM, USA), and data for sites 4 and 5 were 
collected from the meteorological stations of Xinjiang Petroleum Administrative Bureau. At each 
meteorological station, the hourly precipitation was obtained though disaggregation of daily 
observation measured using a standard rain gauge designed by the Chinese Government. 


2.5 Data analysis 


In this study, the differences in BSC variables and soil parameters among different plot locations 
at a macro scale were analyzed using one-way analysis of variance (ANOVA) at the 95% 
significance level. Data analysis was carried out with SPSS version 16.0 software (Chicago, IL, 
USA). The least significant difference (LSD) multiple comparison test was used to determine 
differences among soil parameters and BSC cover values. The data are presented as mean+SE. 
The map of the Gurbantunggut Desert was made by ArcGIS 10.0. The redundancy analysis (RDA) 
bi-plot was created with CANOCO 4.5 and the other figures were created with Origin 8.0 
software (OriginLab, USA). 

CANOCO is one of the most popular programs for multivariate statistical analysis using 
ordination methods in the field of ecology and several related fields (ter Braak, 1988). 
Constrained linear ordination technique RDA was used in CANOCO 4.5 to analyze the influence 
of environmental variables on the cover of BSCs. RDA is a direct gradient analysis or constrained 
ordination method, where the differences in cover measures (response variables) across the plots 
are explained by ordination axes of the environmental variables (Kindt and Coe, 2005). Coarse 
sand (0.500—2.000) and silt (<0.0625 mm) were chosen to represent the distribution of soil 
particle sizes, as they showed the lowest correlation with one another (Lalley et al., 2006a). Soil 
EC and pH were measured simultaneously in the same leach liquor and their values are thus 
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interrelated. Soil pH was significantly correlated with the cover of BSCs in arid environments 
(Lalley et al., 2006a; Li et al., 2010), and then it was chosen to represent the group 
(environmental variables) in our study. It should be noted that in this study, the reduced number of 
environmental variables and the cover of BSCs were standardized before being entered into the 
RDA. 


3 Results 
3.1 Spatial heterogeneity on the distribution of BSCs 


In the Gurbantunggut Desert, there were significant differences not only in the dominant types of 
BSCs but also in the percentage cover at a macro scale. The highest total cover of BSCs occurred 
in the southern site of the desert where moss and lichen crusts covered almost the entire soil 
surface (Table 1). BSCs were slightly less abundant in the central site of the desert and sparsest in 
the eastern site where the covers of cyanobacteria and lichen crusts were 23.83% and 1.54%, 
respectively. In the western site, the interdune surface was covered by 2.16% moss crusts, 28.72% 
cyanobacteria crusts and 59.99% lichen crusts. In contrast, in the northern site, the covers of 
cyanobacteria crusts, lichen crusts and moss crusts were 55.32%, 24.48% and 1.02%, 
respectively. 


Table 1 Spatial characteristics of biological soil crusts (BSCs) in the five study sites 


se Moe Mens Umen Cima anint peo BS 
1 94.1043.17* 18.12+1.89* 73.1842.88" 3.00+0.26° Lichen crusts, moss crusts 
2 89.1743.46* 11.75+1.62° 65.43+42.04*” 8.7541.25° Lichen crusts, moss crusts 
3 80.82+43.38° 1.02+0.084 24.48+1.06° 55.32+2.06* Cyanobacteria crusts, lichen crusts 
4 83.8343.55° 2.16+0.13° 59.99+41.21° 28.72+1.42? Lichen crusts, cyanobacteria crusts 
5 25.17+1.79° 0.00+0.00° 1.54+0.26¢ 23.8341.65° Cyanobacteria crusts 


Note: Sites 1, 2, 3, 4 and 5 are the southern, central, northern, western and eastern sites of the Gurbantunggut Desert, respectively. The 
total cover of BSCs was calculated as the sum of moss crust cover, lichen crust cover and cyanobacteria crust cover. Different lowercase 
letters in the same column indicate significance among study sites at P<0.05 level. Mean+SE. 


3.2 Variation of environmental variables 


Environmental variables (precipitation and soil parameters) differed among the five study sites at 
a macro scale (Table 2). Precipitation generally decreased from the margins to the center of the 
desert (Table 2). The highest mean annual precipitation was recorded in the southern site (133.96 
mm), followed by the western site (109.50 mm) and the northern site (90.27 mm). Moreover, 
precipitation was lower in the central (76.26 mm) and eastern (80.02 mm) sites. 

The spatial variations of soil parameters were not consistent with that of precipitation (Table 2). 
Soil particle size decreased from the northern site to the southern site, and from both the western 
and eastern sites to the central site. The results are consistent with wind directions in the desert. 
Soil pH level and total salt content were highest in the western site, decreased towards the central 
and northern sites, and were lowest in the eastern site. With the exception of the eastern site, there 
were no significant differences in EC, total nitrogen, SOC and CaCO; among the different study 
sites. The highest level of total phosphorus was recorded in the southern site. Moreover, 
significant differences in total potassium were found among the five study sites (P<0.01), with the 
highest value recorded in the central site and the lowest in the eastern site. 


3.3 Relationship between BSCs cover and environmental variables 


The RDA ordination bi-plot of BSCs cover and environmental variables is shown in Figure 2. 
Arrows represent correlations of environmental variables with BSCs cover. The first axis of the 
RDA (eigenvalue=0.56) explained 56.0% of the total variation in the cover of BSCs and 79.7% of 
the relationship between BSCs cover and environmental variables. Soil particle size distribution 
was the most significant environmental variable on Axis 1, represented by coarse sand content 
(canonical coefficient=0.83). The second most influential environmental variable on both Axis 1 
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Table 2 Characteristics of environmental variables (precipitation and soil parameters) in the five study sites 


Site aa ae ~ P Soil pH E a 
1 133.96+14.20° 10.43+2.834 12.44+2.48> 7.99+0.10° 0.09+0.01" 0.39+0.02* 
2 76.26412.378 12.3343.164 14.2743.37° 8.28+0.09° 0.08+0.01" 0.38+0.01° 
3 90.27+12.46° 20.18+3.54° 20.36+3.63* 8.27+0.09° 0.08+0.01" 0.38+0.02° 
4 109.50+13.45° 15.7841.47° 14.77+5.85° 8.43+0.05* 0.09+0.01" 0.41+0.03" 
5 80.02+11.46% 26.39+2.35° 2.85+0.43° 7.90+0.20° 0.07+0.01° 0.30+0.01° 

Site Total nitrogen Total phosphorus Total potassium SOC CaCO; 

(g/kg) (g/kg) (g/kg) (g/kg) (g/kg) 
1 0.21+0.08* 0.34+0.04* 13.65+0.114 1.55+0.59" 12.48+1.01°° 
2 0.14+0.04* 0.21+0.03° 16.55+0.13* 1.21+0.297 12.93+0.28* 
3 0.16+0.04* 0.20+0.05° 15.60+0.32° 1.33+0.30" 12.94+0.56* 
4 0.1740.07* 0.23+0.04°¢ 16.28+0.17° 1.26+0.57° 12.97+40.59* 
5 0.06+0.02° 0.27+0.05° 12.59+0.14° 0.48+0.16° 11.74+0.69° 


Note: Sites 1, 2, 3, 4 and 5 are the southern, central, northern, western and eastern sites of the Gurbantunggut Desert, respectively. EC, 
electrical conductivity; SOC, soil organic carbon. Different lowercase letters in the same column indicate significance among study sites 
at P<0.05 level. Mean+SE. 
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Fig.2 Redundancy analysis (RDA) bi-plot of BSCs (biological soil crusts) cover and environmental variables at 
a macro scale. Moss, moss crust cover; Cyanobacteria, cyanobacteria crust cover; Lichen, lichen crust cover; SOC, 
soil organic carbon. 


and Axis 2 was precipitation (canonical coefficients of —0.68 and 0.58, respectively). The second 
axis of the RDA explained 13.2% of the total variation in the cover of BSCs and an additional 
18.7% of the relationship between BSCs cover and environmental variables. The overall 
ordination explained 87.6% of the total variation in BSCs cover and 100.0% of the relationship 
between BSCs cover and environmental variables. Results shown in Figure 2 demonstrate that 
moss crust cover was significantly and positively correlated with precipitation, but negatively 
correlated with coarse sand content. In contrast, cyanobacteria crust cover was positively 
correlated with coarse sand content. The results suggest that in regions with higher precipitation, 
cyanobacteria crusts would be replaced by moss and lichen crusts. As the content of coarse sand 
increased, lichen crust cover decreased. Lichen and cyanobacteria crusts were also slightly 
affected by soil pH, SOC, soil nitrogen and total salt. Total potassium was shown to be a weak 
explanatory variable with a canonical coefficient of —0.17 on the first axis. 

Results of RDA demonstrate that precipitation and soil particle size (represented by coarse sand 
content) are the main environmental variables that affect the distribution of BSCs in the 
Gurbantunggut Desert. Based on the correlation between the distribution of BSCs and the 
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dominant environmental variables (precipitation and soil particle size), we constructed a 
theoretical diagram to intuitively describe the spatial distribution of BSCs along an environmental 
gradient (Fig. 3). The diagram indicates that the spatial distribution of BSCs is closely related to 
the variations in soil particle size (represented by coarse sand content) and precipitation. 
Cyanobacteria-dominated BSCs occurred in areas with the lowest precipitation where the content 
of coarse sand was highest. Lichen and moss crusts are predominant in the southern site with the 
highest level of precipitation and the lowest level of coarse sand. 
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Fig. 3 Theoretical diagram in describing the distribution of BSCs along the precipitation and coarse sand 
content gradients in the Gurbantunggut Desert 


4 Discussion 


4.1 Impact of environmental variables on the distribution of BSCs 


The impact of precipitation on the distribution of BSCs has been recognized in arid areas 
elsewhere in the world (Eldridge and Tozer, 1997; Belnap et al., 2007; Li, 2012). Zhao et al. 
(2010) suggested that cyanobacteria crusts have higher carbon gains when moisture levels are low, 
for example in sunny, open habitats; in contrast, moss-dominated BSCs have higher carbon gains 
in moist conditions that would be found in shaded habitats. Thus, lower levels of precipitation 
could stimulate the development of cyanobacteria crusts. However, it should be noted that 
cyanobacteria crusts were gradually replaced by moss and lichen crusts as precipitation increased. 
In this study, precipitation was positively correlated with the cover of moss and lichen crusts, and 
negatively correlated with cyanobacteria crust cover. In some areas of the Tengger Desert in Inner 
Mongolia of China, the cover of moss crusts can even reach up to about 100.00% (Li et al., 2010); 
in contrast, in swales (sand dune hollows) of the Gurbantunggut Desert, the cover of moss crusts 
is less than 18.12%, confirming that the relationship between moss cover and precipitation is 
much higher in the Tengger Desert than in the Gurbantunggut Desert. Li (2012) suggested that in 
the Tengger Desert, regions with relatively high precipitation of 200-350 mm are beneficial for 
the formation of moss crusts. Furthermore, lichen-dominated BSCs are usually found in regions 
with annual precipitation of 100-200 mm, and cyanobacteria crusts are restricted to regions where 
the annual precipitation is less than 100 mm. Our study confirmed the results of Li (2012), i.e., 
lichen-dominated BSCs in the Gurbantunggut Desert are found mainly in the southern region 
where the mean annual precipitation is more than 130 mm, and cyanobacteria crusts are restricted 
to regions with a mean annual precipitation of less than 90 mm. However, in the central region of 
the Gurbantunggut Desert, where the mean annual precipitation is about 80 mm, both lichen and 
moss crusts can cover up to 65.43% and 11.75% of the soil surface, respectively. Soil properties 
and dew deposition may contribute to their unexpected successful growth in these regions. 

The RDA results show that soil particle size is significantly correlated with the cover of BSCs. 
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Underlying soils can be modified by the facilitation of dust deposition by BSCs, thus leading to 
an increase in the silt and clay fraction (Li et al., 2010). However, in the Gurbantunggut Desert, 
variation in soil particle size is predominantly due to wind activity at the macro level. The 
proportion of coarse sand is higher and the instability of the soil increases with a corresponding 
increase in wind speed (Wang et al., 2006). As desert soils in the Negev Desert become finer 
textured and more stable, the cover of both lichen and moss crusts increases together with an 
increase in plant species richness (Bowker, 2007; Kidron et al., 2010). In the Tengger Desert, all 
BSCs are distributed in locations where soil particles are of a grain size in the range of 0.01—0.05 
mm, suggesting that particles of this size, albeit smaller than the particle size recorded for silt in 
this study, may provide a suitable substrate and serve as a prerequisite to the formation of soil 
crusts (Duan et al., 2003). In the current study, there was a positive correlation between 
cyanobacteria crust cover and coarse sand content, and a negative correlation between moss crust 
cover and coarse sand content and between lichen crust cover and coarse sand content. 
Cyanobacteria are motile with the ability to move upwards through overlying sediments to reach 
the light levels needed for photosynthesis, while lichens and mosses are incapable of such 
movement and often do not survive being covered by soil (Belnap and Gillette, 1998). In this 
study, few cyanobacteria crusts were found in the eastern site of the Gurbantunggut Desert which 
harbors the coarsest sand. The central site of the desert is characterized by finer textured sand and 
limited precipitation, resulting in high lichen crust cover and moss crust cover. The most abundant 
BSCs developed in the southern site which has the highest precipitation and least coarse sand 
content (Fig. 3). 

In addition to the relationships among the distribution of BSCs, precipitation and soil particle 
size, this study introduces some potential explanations for heterogeneous growth patterns of BSCs 
at a macro scale, such as SOC, soil total nitrogen, total potassium, total phosphorus and soil pH. 
BSCs are carbon- and nitrogen-fixing organisms that strongly influence the chemical 
characteristics of the underlying soil, so the actual impact of soil properties on the distribution of 
BSCs could be greater than the current estimate. A positive relationship among cyanobacteria, 
lichen and soil pH has been documented in numerous studies (e.g., Bowker et al., 2005; Lalley et 
al., 2006a; Li et al., 2010). Unexpectedly, soil pH was negatively correlated with moss crust cover 
in our study (Fig. 2), probably because the impact of soil pH on BSCs depends on the species 
composition of BSCs. Soil total phosphorus and potassium can facilitate the growth of vascular 
plants by increasing shading and limiting access to light, decreasing nutrient availability and 
increasing litter levels, leading to a net negative impact on the growth of mosses and lichens. In 
semi-arid areas of eastern Australia and northeastern Ohio, BSCs cover and lichen diversity 
increase in soils with high CaCO3 content (Grondin and Johansen, 1995; Eldridge and Tozer, 
1997). However, in this study, CaCO3 content had little impact on BSCs cover at a macro scale. 
CaCO; is known to facilitate colonization by cryptogamic plants in arid lands where CaCO3 
content is relatively low. However, only weakly enhanced development of BSCs was recorded in 
our study sites because the soils already contain relatively high levels of CaCOs. 


4.2 Suggestions on the regeneration and rehabilitation of BSCs 


BSCs play many important ecological roles in arid lands, especially obvious for increasing soil 
stability (Eldridge, 1998; Eldridge and Leys, 2003) and nutrient availability (Bowker, 2007). 
Therefore, regeneration and rehabilitation of BSCs are becoming increasingly important for soil 
stability in desert ecosystems. So far, some research has been successfully conducted on the 
artificial culture of moss-dominated BSCs in semi-arid regions (Bai et al., 2003; Tian et al., 2005; 
Rao et al., 2009), and great difficulties in the artificial culture of BSCs were encountered in arid 
regions. However, in the Gurbantunggut Desert, well-developed BSCs have been observed in 
straw checkboards that were fixed for 2—5 years. The straw checkerboard can increase the silt and 
clay content, thereby changing coarse soil texture to fine soil texture (Li et al., 2006; Jia et al., 
2008). Both lichen and moss crusts were developed well in the central region of the 
Gurbantunggut Desert (Tables 1 and 2) where mean annual precipitation is about 80 mm. Soil 
properties and dew deposition may contribute to the unexpected successful growth of lichen and 
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moss crusts in these areas. Regeneration and rehabilitation of BSCs can be addressed in part by 
stabilizing soil and augmenting resource deficiencies such as moisture and nutrients (Bowker, 
2007). Therefore, this research further suggests that if soil is relatively stable and moisture levels 
are appropriate, then the recovery time of BSCs could be hastened. Moreover, the addition of 
phosphorus and potassium has the potential to promote the growth of BSCs and improve the 
overall soil properties. The Gurbantunggut Desert has a huge area of 48x10* km”, so large-scale 
regeneration and rehabilitation programs are not feasible in this way. However, results from this 
study could be incorporated into small-scale regeneration programs along highways and railways, 
and on the margins of towns and villages where the addition of fine grained soil particles (silt), 
moisture and key soil elements including phosphorus and potassium, would be possible. 


5 Conclusions 


The most extensive cover of moss and lichen crusts occurs in the southern site of the 
Gurbantunggut Desert. BSCs are less abundant in the central site of the desert and least in the 
eastern and western sites. The lowest levels of precipitation occur in the central and eastern sites, 
and the highest levels in the southern, western and northern sites. Differences in soil particle size 
appear to be contingent on prevailing wind directions, showing decreasing trend from the 
northern site to the southern site, and from both western and eastern sites to the central site of the 
desert. 

Cyanobacteria-dominated BSCs occurred in regions with lower precipitation and higher 
content of coarse sand. Cyanobacteria and lichen crusts, regardless of whether cyanobacteria or 
lichenes were dominant, gradually appear in regions with higher precipitation and lower content 
of coarse sand. Lichen and moss crusts are predominant in the southern site where precipitation 
level is high and coarse sand content is low. 

Precipitation and soil particle size are the major factors limiting the spatial distribution of BSCs 
at a macro scale. More importantly, the changes in the types of dominant BSCs along 
environmental gradients indicate that increasing precipitation in arid lands has the potential to 
facilitate the growth and development of BSCs. Soil pH and total potassium also explain some of 
the heterogeneous growth patterns of BSCs at a macro scale. The detailed relationships between 
environmental variables and development of BSCs need to be further evaluated by controlled 
experiments in the future research. 
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